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L  INTRODUCTION 


The  solutions  for  space-charge  waves  obtained  by  Hahn  (1939) 
and  Ramo  (1939)  for  a  rectilinear,  drifting  stream  are  quite  well  known. 
Less  well  known  are  the  somewhat  more  mathematically  complicated 
solutions  obtained  when  the  stream  flow  is  allowed  to  be  along  radial 
lines  and  when  the  stream  is  accelerating.  First-order  current  for 
radial  drifting  flow  has  been  obtained  by  Feenberg  (1946).  Current  and 
velocity  in  accelerating  parallel  flow  have  been  obtained  by  Llewellyn 
(1941  ),by  Smullin  (1951),  Tien  (1952),  Hutter  (1954),  Miiller  (1955), 
Higuchi  (1956),  and  Eichenbaum  and  Peter  (1959)*  Accelerating  radial 
flow  was  looked  into  by  Quate  (1952).  The  present  work  attempts  to 
unify  all  of  these  cases  and  a  few  more  in  one  general  solution. 

n.  BASIC  PROBLEM 

As  with  parallel  flow,  it  is  desirable  to  know  how  an  initial 
modulation  is  carried  by  the  stream  from  one  point  to  another.  The 
models  to  be  analyzed  are  shown  in  Fig.  1.  The  mode  with  motion 
along  radius  r  only  and  with  no  variations  of  any  variable  normal  to  r 
will  be  studied.  The  solutions  will  be  for  i^(r,t),  Vj(r,t),  pj(r,t),  and 
E^(r,  t)  for  all  r  and  t. 

in.  GOVERNING  EQUATIONS,  NEW  VARIABLES 

The  linearization  used  in  parallel  flow  is  applicable.  Eulerian 
hydrodynamic al  equations  will  be  used.  Assuming  a  driving  signal  sinu¬ 
soidal  with  time,  the  zero -order  equation  (dc)  and  the  first-order  (ac) 
equations  may  be  separated.  The  dc  potential  V^,  and  hence  the  dc 
velocity  v  ,  is  allowed  to  vary  with  r.  The  velocity  variation  v  (r) 
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will  be  left  arbitrary,  to  be  specified  later  to  suit  specific  problems. 

The  state  of  the  background  is  hidden  in  in  that  the  presence  of 

ions  or  electrodes,  or  just  average  electron  space  charge  is  implied 

by  the  choice  of  v  (r).  Because  the  flow  expands  (or  contracts)  i 
o  o 

and  p  also  depend  on  r. 

°  th 

The  zero-order  equations  are:  the  r  component  of  the 

equations  of  motion, 


dv^(r) 
o  dr 


(1) 


the  definition  of  current 


1  =  p  V  ; 

o  ^o  o 


and  the  continuity  equation,  for  no  variations  normal  to  r, 

V  =  4-  ^  '4' '  “  ■  ''' 

r 

Here  n=0,l,  2  for  planar,  cylindrical  and  spherical  flow  respectively 
The  continuity  equation  may  be  integrated  to  obtain 

i  (T)  =  {7)\(a)  =  {7)"i  .  (4) 

o  r  o  r  oa 

where  a  is  some  conveniently  chosen  radius,  such  as  the  radius  of 
the  inpiit.  Quantities  at  radius  a  will  be  denoted  by  the  additional 
subscript  a.  Similarly 


*  .  oa  .  a  .n 

Po^^)  "  V  (r)  ^  r  ^  ^ 

o 

The  first-order  equation  of  motion  is  just  the  r  component 


equation, 


=  r]E. 


dt  *  1 

The  definition  of  current  is 


‘l  =  •’<.''1  ^ 
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The  continuity  equation  is 

which  is  simply 

1  8  ^**1 


(8) 


(9) 


Several  simplifications  may  be  made  at  this  stage  in  order 
to  reduce  the  algebraic  work*  The  reader  will  recall  that  in  planar 
flow  the  stream  (or  mean)  phase, 


^=w(t-;^)  (10) 

o 

occurs  in  the  solutions  for  i^,  v^,  p^,  and  Thus,  we  anticipate 
that  it  will  also  occur  here,  and  recognize  as 


=  co(t  -  T)  (11) 

where  T  is  the  zero-order  transit  time  from  the  source  to  the  point 
r.  T  is  T(r)  where  changes  with  distance,  as 


T(r) 


■/ 


and 


source 

dT  _ 
dr 


dr 

Vo(r) 


v^(r) 


(12) 


(13) 


Inserting  the  phase  and  a  judicious  choice  on  first-order  qiiantities, 
the  new  variables  are  chosen  to  be: 


i(r.t)  =  i  (r)  +  i  (r,t)  =  i  (r)  +  r""i  (r)  , 

o  1  o  1  -  _ 

v(r,  t)  =  V  (r)  +  v  (r.  t)  =  v  (r)  +  v  (r)e^"  [t-T(r)J  ^ 

O  i  O  1  p  ^ 

P(r,  t)  =  p^(r)  +  Pj(r,  t)  =  p^{r)  +  r ■"pj(r)e^“  [t-T(r)J  ^ 


(14) 

(15) 

(16) 
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(17) 


E(r,  t)  =  E^(r)  +Ej{r.  t)  =  E^(r)  +  r'  ”^j(r)e^“  [t-T(r)]  . 

The  units  of  the  new  variables  are  not  fixed;  for  example*  t  >  has  the 

2  ^ 
units  of  amps/m  for  n=  0*  amps/m  for  n=:l  (to  be  interpreted  as  amps 

per  meter  length/Zir  along  z,  the  cylindrical  axis)  and  amps  for  n  =  2 

(to  be  interpreted  as  the  total  current/ (4Tr)* 

Using  the  new  variables  in  the  first-order  equations,  one 

obtains*  with  exp  jco  (t-T)  understood* 


d^,  dv 

1  .  /V  o  -n;i 
V  — r-  +  V,  — ; —  =  T^r  E, 


(18) 


o  dr  1  dr  “1 

Note  that  the  8Vj/9t  term  was  cancelled  by  part  of  the  (v^  +  Vj^)d(v^+Vj^)/dr 
term*  the  dT/dr  part*  a  saving  afforded  by  insertion  of  the  mean  phase* 
The  remaining  equations  are 


A 

di, 

1  .  0)  A 

57  -}—  ^  +  = 

o 

/>  n  A  A 


(19) 

(20) 


^  A  yi. 

Next*  the  equation  in  i^  and  E^  is  to  be  obtained.  Insert 


from  (20)  into  (19)  and  insert  p  from  (5)  to  obtain 
2  A.  o 

V  di, 

AO  1 

'"l  .  n,  dr 


(21) 


Ja>a  1 


oa 


Differentiate  this  to  obtain 


dv. 


dr 


n 

ja>a  i 


oa 


L  2££' 


dv. 


(22) 


Insert  v  from  (19)  and  from  (22)  into  (18)  to  obtain  the  desired 

^  yK  A  " 


equation  in  i.  and  E 


'  3  ■"'o\ 

A 

,2 

dr 

1 

dr 

111 


oa  *  a  ,n  .  .  _ 

— 7~(  — )  E. 

3  '  r  o  1 


€  V 

o  o 


(23) 
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The  right  hand  side  can  be  simplified  by  recognizing  that  the  ihultiplier 
of  the  displacement  current  is  the  local  plasma  wave  number  (squared) 
as, 


oa 


3 

c  V 

o  o 


2/  V 

2 

-i-  =  Ppi'> 

yo(r) 


(24) 


Hence  (23)  is 

.2a 


d  i,  /  ,  dv  \  di 


(25) 


or,  if  one  likes,  it  is 

dv 

2 


dT^  W  dTjdX-^ 


p(r)(j“‘o%) 


(26) 


The  last  of  the  trio  of  field-driven  equations  to  be  obtained  is 

the  equation  relating  ^  and  E  .  The  steps  in  derivation  are,  roughly: 
2  11. 

take  V  times  the  definition  of  current,  take  d/dr  of  this,  expand  the 

derivative  of  i-v  ,  insert  dL/dr  from  earlier  work  and  i  from  the 
1  o  1  ^1 

definition  of  current?  insert  this  result  into  the  d(^v^)/dr  equation, 

•■n 

factor  out  ^v  and  take  d/dr  of  this,  then  insert  qr  E  for  this 
1  o  1 

derivative.  The  result  is 

dv 


d^(PlVo)  d(?jv^) 


dr 


=  J-E, 

n  1 
r 


dr 


o  I  ^ 

-3—+  1  a 
dr  oa 


V  o 


2  V 

dr  ^^2  J^2  dr  drVVo  ’'o  ‘*7, 


A 

dE, 
,  a  .  n  1 

''‘oa<7' 


(27) 


Equations  (18),  (19),  and  (27)  may  be  used  with  any  form  of 
driving  field. 
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IV.  KINETIC  POWER;  KINETIC  ENERGY;  CONSERVATION  THEOREM 


A  simple  power  theorem  is  readily  derived  for  the  case  of  radial 
flow.  The  equation  of  motion  is 


8v-  av-  av 

1  ,  1  .  o  _ 

“77“  +  V  —  +  v_  -r -  =  Tl  E, 

at  o  ar  1  ar  '  i 


or 


av 

IT  +  ^  <Vo>  =  ’>^1 

ViV  ^  ^ 

Multiply  both  sides  by  (r  and  add  and  subtratt  —  —  (r  i^^)  to 

give 


®''l  .  '”‘l  »<Vo>  .  Vo  3  ,  n.  .  Vo  8  ,  n.  ,  n.  „  „ 

—  IT*  — -57—  ^  —  ‘i’‘~  57  <' 

The  second  and  third  terms  combine;  the  fourth  term  can  be  given  in 
terms  of  dp^/dt  from  the  eqxiation  of  continuity  and  then  combined  with 
the  part  of  the  first  term;  i^  in  the  last  term  should  be  replaced  by 
Jj(t)  -  €^aEj/at  from  the  total  current  equation.  Divide  the  result  by  r*^ 
to  obtain 


+  Ef  +  ^vf+^p,vJ=J,E 


1  2ti  1 


ri 


11 


(1) 


The  terms  here  may  be  recognized  as  the  kinetic  power  density »  the 
electric  energy  density,  kinetic  energy  density,  Wj^,  so  that 

(1)  may  be  written  as. 


(2) 


The  use  of  the  divergence  here  emphasizes  the  meaning  of  conservation, 
that  any  net  power  density  increase  through  a  surface  must  be  balanced  by 
a  decrease  in  energy  density  (for  be  supplied  by  a  source  (Jj?^  0). 
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V.  CIRCUIT  EQUATION;  OPEN -CIRCUIT 


Where  there  is  no  ac  return  path  the  total  current  denisty 

must  vanish,  J, (t)  =  0.  Hence,  the  open-circuit  equation  is 
8E(r,t) 

or 


^^{t)  +  joo  =  0 

Inserting  this  value  of  into  the  current  equation,  (III-25)  produces 

A 

the  homogeneous  equation  in 


2 

+  Pp(r)ii 


0 


(2) 


This  is  a  pleasing  result  from  the  point  of  view  of  the  physics  displayed; 
if  the  acceleration  is  mild,  so  that  the  middle  term  is  relatively  small, 

A 

then  i^  behaves  just  as  in  a  drifting  stream  with  simple  harmonic  motion 

at  the  local  plasma  frequency*  It  is  interesting  to  note  that  this 

equation  is  formally  identical  with  that  obtained  by  Smullin  (1951)  and 

Tien  (1952)  for  the  case  of  parallel  (planar)  flow.  The  effects  of 

2 

radial  flow  are  buried  inside  P  (r)  v/hich  depends  on  n.  It  is  necessary 

z'  P  A 

only  to  solve  the  equation  in  i^  because  v^  and  are  easily  obtained 
from  using  the  earlier  equations.  It  is  worthwhile  to  give  the 
homogeneous  equation  in  velocity^ 


p^(r)  +  — 
p  rv 


dv 
_ o 

dr 


(3) 
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if  only  to  show  that  the  symmetry  among  the  and  pj^  equations 

obtained  in  planar  drifting  flow  no  longer  holds. 


VI.  OPEN  CIRCUIT  SOLUTIONS  FOR  v  (r)  -  r 

o 


m 


The  variation  of  v^(r)  with  r  must  be  given  in  order  to  solve 
(V-2)  or  (V-3).  Among  the  many  possible  variations  of  interest,  a 
common  one  is  that  for  v^  varying  as  some  power  of  the  radius, 


V  (r)  =  (  — )  V 
o'  '  'a  oa 


(1) 


with  corresponding  dc  potential 


vV) 

=  -  -hr 


2t]  'a  oa  'a  oa 


(2) 


These  variations  have  been  used  to  solve  the  planar  flow  model  for 
2m  =  1  (linear  potential  rise)  by  Tien  (1952),  and  more  generally,  by 

II 

Muller  (1955).  The  extension  to  radial  flow  is  relatively  straightforward. 
With  v^  as  given. 


-  dv 

1  _ o  _  m 

V  dr  r  ' 
o 


d^v 


o  _  m(m-l) 
2 


o  dr 

P  pa  r 

(V-2)  and  (V-3)  then  become 
d  X,  ^  di  ^ 

- 1  +  lEl  _JL  +  s2 

’  ^  -  dr  P 


dr 


pa  r 


..a.3m+n  ^  ^ 

(r) 


(3) 
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'^l  ,  2m +  n  ,  fm(m 

j  2  ^  dT 

dr  I  r 


m(m+n-l)  ,  „2  , 

2 -  +  Ppa  < 


a.3m+n|  a 

-)  ]vj=0. 


(■4) 


These  are  Bessel  equations  of  a  standard  from  ^see,  for  example, 
(1945)] 


Jahnke  and  Emde 

l-3m 

«  2  , 
=  r  Z 


A 

Vi  =r 


l-(2m+n) 

2 


and  have  the  solutions 

ib  [±‘''r>'‘]- 

+  (b+l)  J  ’ 


(5) 


(6) 


where  Z  is  a  Bessel  function  of  the  first  or  second  kind,  J  or  N.  The 
quantities  b,  d,  k  are  given  in  terms  of  m  and  n  by 


_  3m -1 

**  2-(3m+n) 

(7) 

_  2-(3m+n) 

=  2 

(8) 

2B  a 

-  _ 

2-{3m+n) 

(9) 

and  of  this  form  satisfy  (111-21)  provided  that  like  signs 
are  chosen  for  both  i^  and  v^.  The  ( +  )  sign  in  the  argument  should  be 
taken  if  the  stream  moves  in  the  positive  r -direction  (expanding  stream 
flow)  and  the  (  -  )  sign  if  the  stream  moves  in  the  negative  r -direction 
(converging  stream  flow)*  For  given  values  of  m  and  n,  the  choice  of 
sign  of  the  order  can  be  made  in  order  to  give  a  positive  value  for  con¬ 
venience. 

It  will  be  convenient  to  use  a  normalized  radius, 

(10) 
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and  a  characteristic  stream  conductance » 


w  2  ' 

7  p  (x)  S  g  X 
•  0)  (x)  ^o'  '  ®a 


(11) 


with  the  dimensions  of  current  density  divided  by  velocity.  The  general  solution 
for  is  a  linear  combination  of  functions  of  the  first  and  second  kind  as. 


and 


AJ^{kx‘^')  +  BN^^(kx^) 


(12) 


kd  -bd+d-1 

— X 


(kx  )  +  BN 


b+1 


(13) 


where  A  and  B  are  arbitrary  constants.  From  (111-21)  and  (1),  (11), 
and  (13)  above,  one  obtains 


(14) 


8a^  L  J 

There  are,  of  course,  a  variety  of  ways  to  specify  boundary  conditions 
in  order  to  fix  A  and  B.  One  way  is  to  specify  the  current  and  velocity 
fluctuation  at  the  initial  plane,  ij(a)  and  v^(a).  This  specification 
gives  the  solution  in  the  form  of  a  transmission  (or  "ABCD")  matrix, 
which  may  be  multiplied  with  similar  matrices,  allowing  a  cascaded 
solution  for  the  more  complicated  dc  potential  distributions.  One 
obtains 


+  j 


[' 


✓ 

^(k)N^(kx‘^)-N^(k)J^(kx‘^)J  a\0  1  , 


(15) 
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Vj(x)a  g(x)  =  -y  X 


_  irk  -(bd+n) 


K 


la 


(16) 


[^b<^>^b+l<^‘^)-Nb<^>'^b+l<^‘^>]  ^“Vla} 


The  bracketed  terms  of  these  eqviations  are  in  the  form  of  the  biradial 
Bessel  functions,  defined  as 


H  (e 

.V 


.?)=|  J 


(e)N^(n-N^  (l)J^(?) 


(17) 


These  and  similar  functions  have  been  defined  and  tabulated  by  various 
authors  and  are  useful  not  only  as  shorthand  notation  but  for  their  pro¬ 
perties  of  differentiation*  integration*  recursion*  etc*  [see  Kino*  (1955 )J . 
Using  the  biradial  functions  in  (15)  and  (16)  and  transforming  back  to  the 
laboratory  coordinate  system*  one  has 


)  kx  ) 


la 


l-m 


) 


,  d-l/2 
kx  exp  j 


^  wa  X  -1 

wt - — - 

V  l-m 

L  oa  J 


(18) 


This  equation  is  most  extraordinary  in  that  (for  v^~r™)  it  covers  planar , 
cylindrical,  and  spherical  streaming,  drifting  or  accelerating  flow 
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VII.  REDUCTION  TO  PLANAR  DRIFTING  FLOW;  n  =  0,  m  =  0 


For  m  =  0,  n  =  0,  (VI-3)  and  (VI-4)  become  the  planar  space-charge 
wave  equations  which  have  solutions  cos  p^r  and  sin  is  instructive 

to  note  that  (VI-18)  also  produces  these  answers.  For  m  =  0  and  n  =  0. 


b  =  -1/2 

d  =  1 

k  =  p  a 

pa 

g(x)  = 

The  biradial  Bessel  functions  now  have  half  integer  orders,  which 
are  particularly  simple,  as 

VT 

=-Ni/2<5)- 

Using  these  relations,  {VI-18)  reduces  to 


cos  6  (r-a) 

pa' 

pa 


j  sin P  (r-a) 
pa 

cos  p  (r-a) 
pa 


(1) 
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Vm.  FASTER  AND  SLOWER  SPACE-CHARGE  WAVE 

AMPLITUDES 


The  first  space -charge  wave  analyses  were  made  by  Hahn  (1939) 
and  Ramo  (1939)  for  a  planar  drifting  stream  and  many  of  their  concepts 
have  formed  a  foundation  for  our  thinking  about  waves  in  electron  streams* 
One  of  these  concepts  is  that  of  ’^faster**  and  **slower'*  space -charge 
waves;  that  is,  the  fluctuation  amplitude  in  space  is  thought  of  as  the 
interference  between  two  constant  amplitude  traveling  waves,  traveling 
faster  and  slower  than  the  average  stream  velocity*  One  would  like  to 
see  whether  this  concept  can  be  extended  to  the  radially  flowing  stream  and 
what  parts  of  the  extension  are  useful* 

The  solution  for  the  planar  drift  space  was  given  in  the  preceding 
section  as 


i^(r,t)=  ^cos 


’  expjj  J^wt-p^(r-a)j| 

Writing  the  cosine  and  sine  in  terms  of  the  exponentials  gives 

+  I  <‘la-''laea'  V'  -jl  } 

ij(r.t)  =  ^ii.+Vj^8,)exp  j  ^t-(P^-p^)(r-a)J 

*l'‘la-''la*a>  “P  >  ["'-(P.+PpaK'-*)] 


(1) 


(2) 

exp  j  wt-p  (r-a)  , 
e 


(3) 


The  first  and  second  terms  of  (3)  are  called  the  faster  and  slower  space - 
charge  waves,  respectively,  with  space -independent  amplitudes  and 
wave-number  s : 


yl  4(iu«u«a>-  P,  =  P.-P 


p» 


(4) 


Is' 


2 '‘la  ’'la«a 


pa 
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The  solution  for  radial  flow  with  v  as  obtained  earlier 

o 

is  repeated  here; 

The  argument  of  the  biradial  functions  is  (k,  kx*^)  and  is  omitted  for 
simplicity.  This  equation  can  be  rearranged  to  have  a  form  a.milar  to 
(2): 


ij(x)  = 


= 


i,  +v,  g  u+jH,  ,)+^{i  -V,  g  -• 

la  la®a  b+1,  b  ^  b,  b  2  la  la®a  b+1,  b 


JH, 


.,b> 


kx 


d-1/2 


(6) 


The  combinations  of  biradial  functions  occurring  in  (6)  become  expo¬ 
nentials  for  large  argument.  For  all  arguments,  these  combinations 
may  be  written  as  exponentials  times  a  varying  magnitude,  as 


(7) 


(8) 


From  the  phase  an  effective  local  plasma  wave  number  may  be  defined, 
for  current,  as 


(9) 


The  bracketed  term  in  (9)  is  unity  for  sufficiently  large  arguments. 

The  faster  and  slower  waves  thus  defined  have  non-constant  amplitude: 
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(10) 


I  ‘la  W I  4'‘la-8a''la'“‘<’‘>-  =  ^a'*'  " 

The  same  method  applied  to  the  velocity  equation  yields  a  different 
effective  local  plasma  wave  number: 

Ck  X  b+l'‘'^b,  b+1  J 


(11) 


=  Pp{x) 


Again  the  bracketed  term  is  unity  for  sufficiently  large  arguments. 
Using  the  kinetic  power  density  as  calculated  earlier, 


_  ^Vo(x)i^(x)v*(x)  J 


(12) 


it  can  be  shown  that 


*"Pj^(x)  =  =  ^kfa  +  ^’kaa 


(13) 


where 


Pj^^{x)  =  Re 


j~Vx)yx)v*(x)j 


(14) 


and  similarly  for  the  slower  wave.  Not  only  is  the  total  kinetic  power 
conserved  in  radial  flow,  but  the  kinetic  power  in  the  faster  and  slower 
waves  defined  by  (10)  is  conserved  separately,  similar  to  the  results  for 
planar  flow. 
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IX.  SOLUTIONS  FOR  3m  +  n  =  2 


Using  3m +  n  =  2  causes  b — >  oo,  d - >  0,  k - >oo,  making 

awkward  the  use  of  the  Bessel  function  solutions.  Such  combinations  of 
m  and  n  include; 


(i)  Planar  space-charge  limited  flow,  m  =  2/3,  n  =  0; 

1/3 

(ii)  Cylindrical  flow  V  ,  m  =  l/3,  n  =  l; 

o 

(iii)  Spherical,  drifting  stream,  m  =  0,  n  =  2. 

Although  some  limiting  process  might  be  used  to  adapt  the  Bessel  function 
solutions  to  fit  these  situations,  a  direct  solution  of  the  current  equation 
(VI-3)  can  be  obtained.  Inserting  2-n  for  3m,  one  obtains, 


d  1, 


dr 


n 


A 

dr 


^2  2 
p  a 


A 

i. 


(1) 


The  solution  to  this  equation  is  of  the  form 

t^(r)  ~  rP 

with  the  determinantal  equation 


+  (1  ~n)p  +  (P„a)^  =  0 
pa 


This  equation  has  two  solutions. 


n-1  ^ 


(2) 


(3) 


n-1  j  n-1  2  ~  ^ 

p=  — ±  V(— '  - 

where  the  quantity  s  has  been  defined  in  an  obvious  way.  Then  one  has 

,n-l  ,  .  .n-1  . 

(o  "i”®) 

i^(r)  =  A'r  +  B'r  (4) 
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A?  (^  +  8)r 


,n-3  . 

r,  1 

+  B'(—  -  S)r 


(^-sP 


(5) 


Vi(r)  =  — ^ 

jwi  a 
oa 

Imposing  the  initial  conditions  at  radius  a  and  transforming  back  to 
laboratory  coordinate  system  gives  the  solution 


V 


g{x)v^(x.,t) 


la 


"  s  .  -s  -  s  -sT  r  L  i  ^  s  -s“ 
X  +x  1-n  X  -X  I  .  fl-n  -  x  -x 

2  2s  2  J  ’^ly  2s  2 


2  ^  S  -ST  r  S  .  -S  ,  8  -8^ 

|l-n  -  X  -X  X  +x  .1-n  x  -x 

^  2s  2  J  |_  2  ■  2s  2 


r, ..  It)! 


(6) 


This  is  obviously  not  an  easy  form  to  use.  Other  forms  can 
be  given  that  are  simpler  and  also  display  the  physical  behavior  more 
clearly. 

The  transmission  matrix  can  be  made  simpler  by  looking  at 

the  nature  of  s.  s  will  be  real  or  imaginary  depending  on  the  magnitude 
2 

of  (6  a)  ,  that  is,  on  the  current, 
pa 

For  planar  flow,  n  =  0,  m  is  2/3,  which  is  fully  space-charge 
limited  accelerated  flow.  This  flow  requires  that 
2 

= I  (’) 


-J  a  A 

a  4 


.3/2 


The  quantity  (P  a)  is  given  by 
pa 


' V’  -  t  h, 

°  via.) 


J  M  a 

a  2  /  a 
a 


(8) 
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Inserting  (7)  into  (8)  leads  to 


(P  a)  =2/9 

pa 


(9) 


so  that 


8=1/6,  p  =  -1/2  +  1/6  =  -2/3,  -1/3 


(10) 


With  s  real,  the  transmission  matrix  elements  A,  B,  C,  D,  are 


A  =  cosh  (s  In  x)  + 


(1-n)  sinh  (s  In  x) 


B  =  jp  a 

pa 


C  =  B  ; 


sinh  (s  In  x) 


D  =  cosh  (s  In  x)  - 


(1-n)  sinh  (s  In  x) 


(lla) 

(llb) 

(llc) 
(lid) 


It  is  seen  that  these  elements  tend  to  increase  or  decrease  monotonically 

with  distance  x  which  is  quite  different  from  the  periodic  elements  of 

oridinary  space -charge  wave  flow* 

For  cylindrical  flow,  n=l,  m  is  1/3  and  the  flow  is  accelerating* 

n-1  is  zero  so  that  s  is  jP  a,  always  imaginary.  With  s  imaginary,  the 

pa 

transmission  matrix  elements  A,  B,  C,  D,  are  using  8=j  ? 


A  =  co.(t  1„  X)  -  ■‘I  ; 


B  =  -jP  a 

pa 


sin(g  In  x) 


(12a) 

(12b) 


C  =  B  ; 


(12c) 
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1  V  Y-  (1-n)  In  x) 

D  =  cos  (?  In  x)  -f  ■  ^  sin  - ^ ^  . 

These  elements,  with  the  trigonometric, functions,  are  periodic 
with  x;  however,  because  of  the  logarithmic  dependence  the  zeros 
of  the  trigonometic  functions  may  be  very  far  apart,  especially  for 
small  1^,  small  currents* 

For  spherical  flow  n-2,  m  is  zero  meaning  that  the  stream  is 
drifting  at  constaht  velocity*  Two  possibilities  arise  here: 

(i)  if  <  0*  5  (small  currents)  then  s  is  real  and  the  mono- 

tonically-increasing  matrix  elements  of  (11)  are  to  be  used; 

(ii)  if  P  a  >  0.  5  (large  currents)  then  s  is  imaginary  and  the 

pa 

periodic  matrix  elements  of  (12)  are  to  be  used* 

The  turning  point  is  at  of  0.  5  which  correpsonds  to  a  critical 

current*  At  less  than  critical  current,  the  behavior  is  similar  to  that 

of  planar  drifting  flow  in  which  space -charge  forces  have  been  neglected; 

the  apparent  reason  is  that  the  flow  expands  sufficiently  rapidly  so  that 

no  debunching  takes  place.  The  critical  current  occurs  where  the  space- 

charge  wavelength,  X.  ,  is  twice  the  circumference  of  the  a  spere* 

pa  — 

At  larger  currents  the  space -charge  forces  have  become  sufficiently 

strong  to  cause  debunching  before  the  stream  has  expanded  appreciably^ 

The  critical  current  can  also  be  associated  with  a  perveance,  as  one  can 
2 

write  (p  a)  in  terms  of  the  total  average  current  sis. 


-I 


pa 


.3/2 


8ir€  ^T) 


(13) 


a  o' 

Using  the  charge -to-mass  ratio  of  electrons,  the  critical  perveance  is 
given  by 
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(amps  /  volts 


(14) 


I 

o 


-33x10 


-6 


a 

(For  example,  this  perveance  corresponds  to  33)i  A  at  1  volt,  or  33 
mA  at  100  volts  or  33A  at  10  kilovolts.  )  For  a  model  in  which  only  a 
sector  of  spherical  flow  is  used,  the  critical  current  and  perveance 
will  be  correspondingly  smaller  by  the  ratio  of  the  area  used  to  that 
of  a  complete  sphere. 

In  the  limit  of  large  radius  of  curvature,  with  small  spacing 
between  input  and  output  and  at  larger  than  critical  perveance,  it  can 
be  shown  that  the  transmission  matrix  elements  are  the  same  as  in 
parallel  flow,  ordinary  space -charge  waves. 


X.  EXAMPLE  OF  CYLINDRICAL  DRIFTING  FLOW  n  =  l,  m  =  0 


The  example  of  a  cylindrical  klystron  is  chosen  to  illustrate 
both  the  physical  and  mathematical  aspects  of  the  rather  involved  solutions 
obtained.  A  sketch  of  such  a  device  is  shown  in  Fig.  2.  For  this  example, 
n=l,  m=0,  which  fixes 

b  =  -1; 

d=  1/2; 

k  =  2p  a 

pa 

g(x)  = 

For  the  sake  of  convenience,  as  there  is  no  critical  current  or  perveance, 
P^a  may  be  chosen  arbitrarily;  one  easy  choice  is  to  make  =  1/2  and 

thus,  k=l.  The  relation  between  input  and  output  quantities  from 
(VI-18)  is. 
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Figure 


g(x)Vj(x) 


la 


exp 


\  Via 


[-j(x-l)] 


(1) 


(Use  has  been  made  of  the  recursion  relations  among  the  H  functions.  ) 

.V 

For  velocity  modulation  only  at  the  a  (initial)  plane 


Vj(x) 

'"la 


[M 


f 


=  jHi  od.y^D 


(2) 


and 


\(x) 

vIT 


^1, 1^^'  ’ 


(3) 


The  pertinent  biradial  functions  are  shown  in  Fig.  3  as  a  function  of 
X.  Note  that  the  velocity  and  current  density  are  periodic  as  with  planar 
drifting  flow,  but  here  the  velocity  maxima  and  current  minima  are  not  quite 
aligned,  except  at  large  distances.  The  bunching  and  debunching  is  much  the 
same  as  in  ordinary  parallel-flow  klystrons.  However,  the  total  ac  current 
maxima  will  not  stay  constant  but  will  grow  with  the  distance  as  can  be 
shown  using  the  asymptotic  formula  for  For  ?,  5 - ^  od 


^  sin  (?-?)  - 

~7f? - ^ 

Thus,  for  large  values  of  k  and  kyx,  the  current  (~  i^(x)  x)  is. 


(4) 


^l(x)  X  ^  sin  (k/x  -  k) 


Via 


f7 


(5) 


3/4 


so  that  the  current  maxima  growth  as  x  for  expanding  flow. 
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(1,/x}  ^velocity 
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For  current  modulation  only  at  the  a  (initial)  plane 


(v.  =  0),  one  has 

la  ' 


,  ,  v,(x)  g 

IT- 

la  ’  la 


(6) 


and 


i,{x) 

T-  =  '“o,!"’ 

la 


(7) 


The  pertinent  biradial  functions  are  shown  in  Fig.  4.  The  same  comments 
apply  as  in  the  case  of  velocity  modulation. 

An  interesting  singularity  in  the  charge  density  occurs  under 
certain  conditions.  Recall  from  planar  drifting  flow,  with  only  the 
faster  wave  excited,  that  the  ac  charge  density  vanishes  when  excited 
at  w  =u>p.  In  the  cylindrical  drifting  flow  with  the  stream  expanding  , 
the  local  plasma  frequency  decreases  with  increasing  x.  Thus,  by  modu¬ 
lating  at  r  =  a  with  a  frequency  c*>  <  and  letting  the  flow  expand, 
there  will  be  a  radius  at  which  co  =w^(r).  The  question  is  whether  the 
ac  charge  density  can  be  made  to  vanish  at  this  radius.  For  the  model  at 
hand,  with  n  =  1,  m  =  0,  the  charge  density  is  given  by 


Pi(x) 

lj(x) 

Vj(x) 

p 

i  v 

^oa 

oa 

oa 

Suppose  that  at 

some  radius  x 

xa 


(8) 


il(Xc)  =  Vj(x^)g(x^)  (9) 

which  says  that  the  velocity  and  current  are  in  phase  at  this  point,  and 
implies  the  existence  of  a  **faster  wave**  only  at  this  points  Using  this 
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I 


-26 


Figure  4 


value  for  the  velocity,  the  charge  density  becomes,  at  x^, 


\(x^) 


(1  - 


..  P 


oa 


oa 


0) 


(10) 


Hence,  the  charge  density  can  be  made  to  vanish  at  the  point  by 
choosing 


w  =  w  (x  )x 
p  c  c 


(11) 


This  zero  in  charge  density  comes  about  not  by  wave  interference  but 

by  a  local  physical  requirement  of  zero  charge  density*  In  order  for 

the  modulation  at  x  to  be  that  specified  by  (9)»  the  modulation  at  the  a 
c  ^ 

(initial)  plane  must  be 

■-H,  (i.y^  +  jH  jd.  f^)n 

- E - hi - £_  (12) 

L-,  o,  I  c  o,  o  c  ^ 

The  bracketed  term  is  unity  both  for  x  =  l,  and  large  values  of  x  and 
tends  to  have  magnitude  near  unity  for  values  in  between;  thus,  requiring 
the  ''faster  wave"  only  at  x  =  implies  that  something  like  the  "faster  wave" 
exists  for  all  other  x* 


^la  =  ^la^la 
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T.*}.*  irili  Al  1'  I  itiiii'rnia 

J  tit  'h  V  Hi  --'.ir'  li  Cl  filler  (NASA) 

AT  TN  IiM  ni  ii.il  I  ihrary 
l.ar*i.1<  \  Al.'  C.ilil'ii  .iia 

L.vti.  . .  I.  rm.i-r  (NASA) 

ATIN  I . .  »'  LiCr.iry 

2l<K)ri  III  iMikihirV  Rimi'I 
Clt*vr-litn«i  t>*  I'l 

Tiiiflilirrl  h|if  ■  Flijilii  Ci-nitr  (NASA) 

A  I  TN  )>'<  . . .  Library 

ij rrrniit.li .  Mn rvl.m'l 

Oi-.if'it-  C  R'.n-h.ill  ‘.jwti  i‘  Flijihi  Canter  (NASA) 

AT  )  N  I . .  l.'i  I.  I  V 

Hi'd.l'itiv  Ar  .iMi.il.  /'i.i'h.rna 


a.,ll  „i-.  M,iii-.n  (NAs.'.l 
AT  TN,  li-  111*'*  al  I  iiirtry 
WUlop-*  Itland.  Vir|{ini,i 

Inaiiliiii-  ot  At  rir.aulicAl  Sciancaa 

2  Eaal  b4(h  Slrett 

Nan  Yorli  21  New  York 

Ap|.tied  Wchanita  Review* 
Souihweel  Rcreanh  Inalitiite 
2100  Culehra  Roail 
San  Anl  inm  b.  Ti'ea-- 

l.iiida  Hall  Library 
ATTN'  tWiimenl  Dtvitmn 
ciOi  Chtrry  Sirori 
Kan.ii*  Ciiy  10,  Mi*  .nun 

AFOSR  (SRAM 

ATTN  Tichiiicai  Library 

llotloman  AFO.  New  Meeicu 

AFSWC  (SWOl) 

KirilandAFA.  Ne*.  Meaieo 

Advam  rd  Reiean  h  ProjecU  Agency 
Waihin^l'Hi  2S  D.  C 

Rnnd  Corporal iim 
1700  Main  Siren 
SanU  Mcniva,  California 


(UncUavtli 
Report  e) 


(ttrc  laneificd 
Rep<  rli) 


(ItncUetiRrd 

Report*) 


Cliairman 

Canadian  .*i>lPt  Staff  (DRF/DRISI 
2410  Ma**ayliwaeit«  Aveawa,  N.  W. 
Watbinptr.r  2i|  O.  C. 

Otfiver  in  CtarR* 

OtiKv  uf  Naval  Re«earch 
Nav,  N».  100 
Fleet  P*  »l  0*f*c.* 

Ne*  Yr  rll  N.w  Y-rh 

OKiif  .»1  lethrm  •!  Vrvue* 
Dwpaetniant  of  C>  mmerre 
Tachniral  Report*  Brant  h 
ttarlitRii'ii  21.  £  C 

tA  D  .Vdtmb.  Heal 
Ament  III  Syateme  Inr e r )■  irate d 
l#.21  t  I2b'»|.***i 
Hawthar.ic,  Caltl'.mta 

Or,  Harold  Ciarer 
Oflite  'tf  Naval  Rraaarrti 
Wantimptor  21.  0.  C 

Hrolo**u'  Clw'lei  Townaa 
Oepannaei.t  •>(  »'hvei*a 
f  eliimnia  •JinvcreHy 
Ntw  Turk  2  7  Ne«  Vorb 

l*rntr«>t>ir  Harvey  brixRia 
Oepartmmt  of  Pliyetc* 

Harvard  Daiveraliy 
Cainbndpe  12,  MaaaaclMaatt* 

Rrriaaai.r  P  Koai  li 
naparinwnt  <•*  Phyaica 
Ctlwmbla  tinivrraily 
New  Vorb  27.  Nee  Ynrb 

Prefaaaor  H.  ttlieirberscr 
Hepartment  of  Phyatra 
Harvard  Uaivrraity 
CaiaSrtdte  M.  MaraachaaMt* 

Or.  IrviM  Rime 
Olftca  at  Naval  Reeeapclt 
MS  Braadway 
Na»  Tara.  Naw  Torli 

W.  A.  lUBaMplIk.  Miner 
LarbReei  Aimall  Carperaili 
Rheattae  aad  Rpia  DinaH 
Tackalral  MrataHli  Caiar 
1291  Maaeeer  Mraat 
PaloAlu.  Calllaria 

Mr.  B.  Ohraaa 
•Mrry  OiretaMe  Caaipaay 
TaSa  MvteMa 
MeU  BiaiMn  t  R4« 

Great  Nrrit.  Ne*  Verb 

OOFL  (OROTL  Ml) 

WaetMnglti*  29.  D.  C. 


(tin'  lareilied 

Raparta) 


Oral.  J  N.  MidlUaa.  ir. 

CRMnni.  OapRrtiaat  ef  Bactrteal  Batnerfag 


ADDRESSEE 


NO  COPIIn 


The  Ohio  State  Untveraify  Antenna  Lahoraiury 
2024  Neil  Atenue 
C'llomhaa  10,  Ohm 
ATTN:  Setgrily  Oliiter 

Miiliir  (-.eneral  Cateni.r.i  M;«filrnf2rii  Filhj 
Crniral  Tift  hniio  da  Aeronaullta  (CTA) 

Sio  Jo(^  d’^e  Campu* 

SCo  Paulii,  r.raeil 


Hiirdiiv  Univereily 
Lalavnie,  Indiana 
ATTN;  RithardL.  Fonkho-i^er 
Enaineerina  Librartan 

Pri>(eii*Ar  J.  Van  Bladrl 
Fleclncal  Enptnerrlnp  Department 
C'lllear  t^f  F.n2inerriii2 
Unitersily  ol  Wi*«on*in 
htadiBon  S.  Wtecaiein 


pr't(e«k>e  R.  l.«iw 
Duperinieni  <il  Fhync* 

The  Hebrew  linlvervny  .»f  Jcruealrm 
Jrruealem,  l*reel 

Univerrity  of  (I  I  mole 

Department  ol  ClrcirKai  Enalnvenna 

Urbane,  lllinoi* 

ATTN.  H.  V.m  Forrrter 

Ohio  Sttle  Univeriitv 

Deperiment  >il  EUctrical  Eniiirmny 

ColuriibttB.  Ohio 

The  UniverBity  u(  Michigan 
Department  of  CIrttrtcal  Enrlncrrina 
Eli'ilron  Phyeic*  Laboratory 
Ann  Arbor,  Mxhigrn 
ATTN:  Prolreeiir  2.  E.  Rowe 


Stanford  Inlvereity 
Electronic*  Rcaearrn  Laboratory 
Sianf'ird.  California 

Marearhitaeite  Inriilule  of  Teihnolngy 
Reeearch  Labe,  of  Clectroatci 
Room  20R-22i.  Doc ument  Office 
Cambridge  19,  Maatachueelt* 

ATTN:  J.  H.  Hewitt 

Harvard  l'nivcr«it^ 

Cruft  Lmboratary 
Cambridge  11.  Mae*aihu»ett* 

ATTN;  Tethnicil  Report*  CoUecIton 

TvchnKal  Information  Ltbrartrt 
Rell  Trlvphone  Laboraioriea.  Inc. 
Wbippany  Laboratory 

Wltippany-  Now  Jereey 

ATTN;  Trrhnutl  Rrpi.rt*  LibrariAii 

Caltfornta  Inetituta  of  1  .-chnoloiv 
Paeadrna  -I,  California 
ATTN;  R.  Could 

.Sylvanta  Claitnt  Company 
btountain  View.  CaHfnrnia 
ATTN:  O.  K.  CoodfMn 

PrUteaor  Willtam  H.  Surder 
F.>ectr<tal  Enftnrertay  IhrpartmH 
Priacrtmi  UatveraWy 
PrlMoitm.  Naw  Jtreay 


Goaoral  BIrci rU  Company 
Electriiii  Tebe  DIv.  of  the 
Tho  Rpnll* 

Sclwaevtadl,  Hew  Yora 
ATTN:  E  O  McArtbar 


Reeearch  Lah. 


Ceteml  llecirtt  Company 
bBaiila  ana  Sp**^*  Vohicia  OryatSKioM 
It'll  Chaafavf  Preet 
PhllaUphia  4,  Hrnaiy’eanta 
ATTN:  Arrodynamite  EngNtaerlAji  Operatto* 
M  F.  Pach. - 


Hwfhra  Atrrraft  Company 
noreaca  at  Toala  Sirrof 

Calver  C«y.  California 

ATTN;  Oacaimnta  flraap  Bid*.  S,  Rm  K«9 


RCA  Laboraionae 
Prtfcotoa.  Naw  ieny 
ATTN.  Or.  «  M.  Wehoiar.  Oirertor 
Eloftrontca  RoaMrch  Laba. 


Vartta  Aaeoctatoa 
SU  MaaM  Way 
PUa  AHa.  Callfaeau 
ATTNt  Taelaiical  Lthrary 


t  nactrw  CorporatUia 


Battranie  Taba  Obvtaii 
P.  a  Bl  AM 
Bltnfra.  Mir  Terh 
ATTN;  Sir.  r  ■* 


Prafaaaoe  Ealirab  Kapawtal 
OnlvaraMy  al  liaffcata  Callfaiia 
•elMwl  af  EngNMaHae 
Dipimai  al  BactrlaM  BatMaMMB 
OMearnRy  MA 
Loa  Angalaa  9,  CaltfoenM 


piilaia  VirlTT  Reaaareh  LaborMortee 

OelMaM.  Cnllfarnta  I 


HaadgaarMra 

Air  Pnraa  OIfka  al  BcteMiflc  fUaadrtfc 
Offtea  M  AaeaafMe  Maaanh 


_ _ ^_J|29.  ».  C. 

ATTWt  MspakallC.  MaatMipi 
Oamipl  PRyalca  DMatan 


ar 


